Quantum Field Theory and the Internal States of Elementary Particles 
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A new application of quantum field theory is developed that gives a description of the internal 
dynamics of dressed elementary particles and predicts their masses. The fermionic and bosonic 
quantum fields are treated as interdependent fields satisfying coupled quantum field equations, all 
expressed at the same space-time coordinate. Quantization is realized by expanding the quantum 
fields in terms of fermionic creation and annihilation operators. This approach is applied in a QCD 
description of the light quarks with a zero Higgs field. Originally massless and pointlike, an isolated 
quark (described in its own center-of-mass) acquires mass and a finite extent when treated as an 
interacting system of quark and gluon fields. The binding mechanism of this localized system has a 
topological character, being a consequence of the non-linear nature of QCD, while being insensitive 
to the magnitude of the coupling constant to lowest order. To prevent this system from collapsing 
general relativity is introduced. The quark stabilizes at a radius of 8.8 Planck lengths and acquires 
a mass of 3.2 MeV, in remarkable agreement with accepted phenomenological values. It is suggested 
that the two higher generations of quarks are associated with the other two real solutions of the 
Higgs field equations. 
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I. INTRODUCTION 

The standard model has been extremely successful 
in describing scattering phenomena between elementary 
particles. Nonetheless, it is clear that this picture of Na- 
ture cannot be complete. Most importantly, the standard 
model contains too many free parameters to be called 
fundamental, and can thus best be seen as an effective 
theory. In modern treatments the masses of the standard 
model fermions are supposed to arise from the couplings 
of quarks to various Higgs fields, however, so far this has 
not lead to a reduction of the number of parameters. Of 
the 19 free parameters of the standard model, 13 belong 
to the Yukawa/Higgs sector [l|. These parameters are 
put in by hand, and there is no explanation for the huge 
hierarchy of masses ranging over 6 orders of magnitude 
(not including the neutrinos) . Further extensions beyond 
the standard model often lead to a greater - rather than 
a smaller - set of unknown parameters. 

Given this situation one may wonder whether it is not 
possible to underlie the standard model by a theory at 
a more fundamental level that evolves into the standard 
model through the construction of effective solutions of 
the bare degrees of freedom. Such an underlying the- 
ory should contain many of the elements of the standard 
model but have a simpler structure, i.e. it should embody 
a certain unification of the diverse elements of the effec- 
tive theory. The simplicity should imply fewer degrees of 
freedom and fewer parameters at the fundamental level. 
The new application of QFT methods, introduced in this 
paper, enables such a link between fundamental and ef- 
fective levels by its ability to construct massive dressed 
elementary particles from more fundamental bare mass- 
less particles. 

An important element of such a fundamental theory 
could be the reduction of the three generations in the ef- 
fective theory to a single one in the fundamental theory. 



Since the three generations in the standard model are 
not associated with any additional symmetry or interac- 
tion, three generations may well be the result of multi- 
ple solutions at a more fundamental level. This possible 
scenario would gain in standing if we could explain the 
emergence of (exactly) three generations from the exact 
level. We suggest that a single fundamental Higgs field 
might explain the threefold split. The Higgs field equa- 
tions would typically have three real solutions, whose 
classical approximations could be denoted by <fi = and 
<p = ±^/|/i 2 /A|. In the quark-Higgs dynamics the pos- 
itive and negative solutions lead to distinct solutions, 
so that these three solutions would indeed lead to three 
distinct fermion sectors. Since the nature of the Higgs 
Lagrangian at the fundamental level is not (yet) unam- 
biguously established and a non-zero Higgs field would 
complicate our basic model considerably, we limit our- 
selves in this paper to the trivial Higgs solution, namely 
<f> = 0. However, this also eliminates the Higgs parame- 
ter n from the model, so that it is unclear how this basic 
theory acquires a scale. The question of the basic scale 
parameters in Nature has been considered in the context 
of cosmology elsewhere 0. We find that for the light 
quarks general relativity has to be introduced to ensure 
the existence of the quarks, so that for this basic gen- 
eration the cosmological parameters reviewed in Q also 
play a role. The presence of different scale parameters at 
this fundamental level could give a possible explanation 
of the hierarchy problem in physics, as it would allow for 
the existence of a range of different scales at a very basic 
level. 

Which further ingredients would survive in a funda- 
mental theory underlying the standard model? It is to 
be expected that the SU(3) theory of strong interactions 
will survive, as this theory is unbroken at the standard 
model level and is not characterized by a complex mul- 
tiplet structure. The complexities of the multiplet struc- 
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ture and broken symmetries in the standard model origi- 
nate from the electro- weak interaction and it is not clear 
which of these complexities emerge in the transition from 
fundamental to effective theory, or which ones are al- 
ready present at the most fundamental level. Therefore 
it is natural to limit oneself originally to QCD. Clearly, 
so long as one does not consider the electro-weak inter- 
actions one cannot describe leptons as effective systems, 
so we also focus on quarks in this paper. Being a non- 
linear theory, one expects that QCD can supply the nec- 
essary self-consistent binding mechanism for the effective 
quark system, recalling the important role non-linear dy- 
namics plays in soliton particle models in hadron physics 
0. By adding the electromagnetic interaction at a later 
stage one can break the isospin symmetry of QCD and 
try to explain the mass differences between up and down 
quarks. Being a linear theory, the electromagnetic force 
can be treated perturbatively and is not expected to in- 
terfere with the expected non-linear binding mechanism, 
so that it can be ignored initially. 

As stated before, the possibility to link the dressed par- 
ticles and their internal properties and masses to more 
fundamental bare particles relies on a new methodol- 
ogy which emphasizes the use of the field equations in 
QFT. The standard applications of quantum field the- 
ory (QFT) deal with scattering processes expressed in 
term of propagators and vertices, yielding the well-known 
Feynman diagrams. The dressing of fermions is then ex- 
pressed in a series of time-ordered interaction diagrams 
modifying the bare propagator. Most of these diagrams 
are divergent, which has necessitated the renormaliza- 
tion methods to arrive at meaningful physical results, 
although outside inputs are still required to fix the phys- 
ical parameters. This dressing process does not yield any 
internal properties of the particles, due to the fact that 
the fermions interact with external bosonic fields. The 
dressing process which is introduced in this paper is of a 
very different nature as the fermionic and bosonic fields 
are treated as mutually dependent fields, the dependence 
being controlled by the coupled field equations. Hence, 
the boson fields are not treated as external fields and all 
fields and interactions refer to the same space-time coor- 
dinate. It is this treatment that enables one to describe 
the internal properties of elementary particles. 

It is clear that the field equations play a dominant role 
in the new formulation. This is not unlike most other ar- 
eas of physics. However, in QFT their use has been very 
limited to date. Often the field equations in QFT are de- 
noted as classical, and indeed they are often used as such. 
For example, the scalar Higgs field 4> satisfies a classical 
field equation with various (constant) classical solutions. 
The "full" solution is then expanded around this classi- 
cal solution, and only the additional fields are quantized. 
Despite their classical nature the constant solutions are 
attributed physical meaning, and their value is often de- 
noted as the vacuum expectation value of the Higgs field 
[l|. This is a dubious use of the term expectation value 
as it suggests that one takes the expectation value of an 



operator, which one does not. It is not sufficient to call 
something an operator, one also has to treat it as such. 
The expansion around classical solutions is also common 
in the soliton and instanton problem Q. Here the basic 
solution also satisfies an approximate non-linear classical 
equation of motion, and quantization is imposed by lin- 
earizing the remainder. We feel that these procedures are 
fundamentally flawed. If classical solutions have any role 
to play in QFT then they should emerge as an approxi- 
mation to the quantum operators, not lie at the basis of 
the quantization procedure. Quantization should be done 
up front, not as an afterthought. Hence, we promote a 
much more rigorous quantum-mechanical use of the field 
equations. Before solving the field equations one should 
turn them into quantum equations and ensure that the 
solutions are quantum field operators. 

A rigorous quantization procedure can be realized by 
turning the classical fields into quantum fields by expand- 
ing them into creation and annihilation operators, the 
latter satisfying well-defined (anti-) commutation rules. 
The insertion of quantum operators in the field equations 
also means that the order of the fields in the field equa- 
tions can make an important difference, a phenomenon 
well-known in most quantization procedures. Hence, be- 
fore solving the field equations one must establish the 
correct order of the fields. This is usually possible by de- 
manding consistency, by imposing symmetries or through 
other physical requirements. At this point it is of inter- 
est to review the standard quantization procedure for a 
linear field theory such as QED, since this may eluci- 
date why classical solutions have played such a promi- 
nent role in traditional quantization procedures, a role 
which should be avoided when one deals with non-linear 
equations. To derive the electromagnetic field operator 
in QED one first determines the classical solutions of the 
classical linear equation (plane waves in the free case). 
One then quantizes the field by appending creation and 
annihilation operators to the classical solutions. This 
well-known procedure leads to the correct quantized elec- 
tromagnetic field, and therefore may well be responsible 
for the practice to use classical solutions in the quanti- 
zation procedure. However, we would have reached the 
same result by inverting the order, namely, by first ex- 
panding the field in terms of creation and annihilation op- 
erators of (as yet) unknown functions, and then solve the 
(quantum) equations of motion. This procedure would 
yield the same plane wave states. However, by choosing 
this latter order one guarantees up front that the solu- 
tions are quantum solutions. For non-linear field equa- 
tions it makes a huge difference which order is applied. If 
one solves the non-linear equations classically, and then 
appends creation and annihilation operators to these so- 
lutions, then the resulting operator is generally no longer 
a solution of the quantum field equations, as the oper- 
ators do not behave like c-numbers. Hence, one must 
first introduce the field operators and then construct op- 
erator solutions to the field equations. If the resulting 
profile functions, which multiply the creation or annihi- 



3 



lation operators, are identical - or close - to the classical 
solutions, then this could justify a classical approxima- 
tion. However, such a situation only arises if the resulting 
field operator acts like a unit operator in operator space. 
In general this situation docs not apply and classical so- 
lutions cannot be expected to play a fundamental role in 
non-linear field theories. 

As indicated above, creation and annihilation opera- 
tors play a central role in the quantization procedure. 
The operator nature of the field now even becomes part 
of the solution process and in general the simple linear 
expansion of the fermion and boson fields in terms of cre- 
ation and annihilation operators is no longer a solution 
of the full set of field equations. One can start with the 
simple linear expansion of the fermion field, however, af- 
ter further iterations a much more complex structure will 
evolve. For example, to lowest order the gluon operator 
becomes a bilinear operator in terms of the fermionic cre- 
ation and annihilation operators. This bilinear expansion 
of the gluon field in terms of fermionic operators, instead 
of in single bosonic creation and annihilation operators, 
does not imply that gluon fields have lost their identity. 
In fact a similar expansion of the gluon field (or the pho- 
ton field for that matter) in terms of quark operators is 
possible when it is considered an external field (and thus 
fully independent field) and satisfies the homogeneous 
gluon field equations. So such an expansion in fermionic 
degrees of freedom is perfectly consistent with the idea of 
independent gluons, soft gluons and glueballs. However, 
in the self-consistent case, which we are considering here, 
the gluon and quark fields are interdependent and neither 
of these fields can be considered as independent, although 
it still makes sense to talk about the fermion and gluon 
degrees of freedom. One should also not interpret this 
formal expansion of the gluon field as a composite model 
of the gluon. It is known that the quark- antiquark in- 
teraction is repulsive in the color octet gluon channel, 
however, forces play no role in the operator expansion 
and this fact plays no role in the expansion. 

By further iteration one can construct the whole op- 
erator structure of the exact quantum fields. Ultimately 
one arrives at an infinite sum of operator terms, each 
additional term containing a higher power of creation 
and annihilation operators. Despite the elaborate op- 
erator structure of the fields and the equations it is pos- 
sible to construct a closed operator solution. In Section 
IITI and Appendix [A] this complete operator structure of 
the fields is derived and is expressed in an elegant con- 
cise form, involving a vacuum projection operator. By 
factoring out the quantum operators the quantum field 
equations can then be reduced to ordinary coupled equa- 
tions, similar in structure - but not in content - to the 
classical equations. The elegance of this reduction is fur- 
ther underlined by the fact that the number of profile 
functions remains finite, despite the infinite nature of the 
expansion. The dependence on spin, color and isospin 
can also be factored out, leaving in the end a set of cou- 
pled differential equations for c-number profile functions. 



The number of profile functions depends on the symme- 
tries present: for our QCD formulation there are only 
5 gluon profile functions. But if symmetries are broken 
(like isospin in QED), more profile functions are needed. 
One of the great unexpected benefits of this quantum 
reduction is that the resulting non-linear equations are 
more amenable to analytic solutions than the original 
classical field equations. 

After further extensive manipulations of the scalar cou- 
pled equations, two simple order 0(a®) solutions emerge. 
The first one is an uninteresting "trivial" solution, in 
which all effective potentials vanish. The second one fea- 
tures non-zero potentials that become singular at a fi- 
nite radius tq. This absolute confining mechanism of the 
quark within its self-generated bag is of great beauty and 
is made possible by a unique QFT mechanism not seen 
before in particle theory. Since this solution is enabled 
by the non- linear nature of the field equations and is in- 
dependent of the coupling constant, it will be denoted as 
a topological binding mechanism. However, this solution 
still features some other problems. First, in the absence 
of a Higgs field (we set <fi — 0) this QCD model does 
not contain a scale, so there is no indication whether this 
model refers to a real particle. A second problem is that 
the expectation value of the energy operator is negative, 
so that this energy cannot be identified with the mass 
of the quark. In addition, a system with negative en- 
ergy will collapse to a point as this minimizes its energy. 
These problems have an amazing solution. By linking 
QCD to general relativity (GR), we can exploit the fact 
that GR contains a definite scale and that gravity will 
act as a repulsive force if the energy of the system is 
negative. Hence, after the system has contracted to the 
Planck scale the repulsive gravitational force becomes big 
enough to halt the collapse. By adding vacuum energy 
to the system we can compensate for the negative internal 
energy, and eventually obtain a net positive energy that 
can be identified with the mass of the system. Amaz- 
ingly, this mass comes out to be slightly over 3 MeV, 
which is close to the phenomenological values found for 
the light quark masses. Hence, this solution can be iden- 
tified with the light quark doublet. This solution also 
resolves the apparent conflict between the smallness of 
the quark mass and its pointlike nature. Normally, a 
mass of 3 MeV would correspond to a ridiculous size of 
65 fm, but our result of 8.8 Planck lengths lies well within 
observable limits. 

Now that we have explained the basic elements in our 
formulation of the internal bound state dynamics of QFT, 
it is appropriate to mention earlier approaches in QFT 
that tried to describe localized states. Soon after the de- 
velopment of the Dirac equation, this equation was used 
to carry out bound-state calculations, for example by de- 
riving relativistic corrections to the non-relativistic hy- 
drogen atom calculations. The success of these calcula- 
tions already suggested that QFT can be used fruitfully 
in bound-state problems, however, the nature of these 
Dirac calculations (the proton is not treated as a quan- 
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turn field) shows that such a calculation cannot really 
be called a genuine QFT calculation. Historically, the 
electron also has been subjected to bound-state calcu- 
lations of quantum nature. Schweber Q reviews efforts 
to describe the electron as a bound state in the period 
1940-1955. Dirac i(| made an attempt to describe an 
electron as a particle with a finite, perfectly conducting, 
surface in 1962. However, these efforts were hampered 
by the fact that the electroweak and strong gauge theo- 
ries had not yet been developed at the time, while one 
can expect that the non-linearity of these non-Abelean 
gauge theories is essential for the construction of bound- 
state solutions. Other methods have been developed to 
handle bound-state problems within standard QFT, for 
example, the Bethe-Salpeter equations (El, 01, [Icj . 
[IH). the Blankenbecler-Sugar equations and the 

Dyson-Schwinger equations [13|. However, none of these 
theories could describe the internal properties of the el- 
ementary particles as the bosonic fields were treated as 
external fields. More recently, lattice gauge calculations 
have been applied to nucleons considered as assemblies 
of quarks However, these calculations are also not 
designed to give insights in the quarks and leptons them- 
selves. 

Historically, our efforts to solve the non-linear QCD 
field equations self-consistently started out as an attempt 
to describe the binding of three quarks in a bag, expect- 
ing to develop an alternative to the bag models popular 
in the eighties fl5j . Solutions were found fairly quickly 
[l6j . however, their physical interpretation left much to 
be desired as the mass/energy scale of the light quarks (a 
few MeV) corresponds to very extensive objects, which 
are physically unacceptable. Also, only particle states 
were considered initially, as was common in the field of 
intermediate energy physics [l7j and other papers deal- 
ing with pion corrections to quark bag models [l5j |. It 
soon became clear that it is difficult to apply such an ap- 
proach without anti-particles consistently, and this initial 
approach was replaced by a formulation in which all fields 
are expanded in particle and anti-particle creation and 
annihilation operators. When this formulation lead to 
the exact operator solution it also became clear that this 
approach cannot be used for composite systems such as 
three quark bags, as the operator solutions only survive 
if they operate on one-body state vectors. Hence, the 
attention switched from describing composite systems to 
elementary particles and this led to the realization that 
these new QFT methods can be used to give a descrip- 
tion of the internal dynamics of elementary particles, in 
particular quarks. 

The outline of the paper is as follows. In Section |TT] 
we discuss the equations of motion and derive the form 
of the field operators in terms of quark creation and an- 
nihilation operators. The formal solution of these op- 
erator equations is presented, allowing the reduction of 
the field equations to ordinary differential equations. Af- 
ter defining the c-number wave function components and 
c-number gluon profile functions in Section IIII1 we de- 



scribe the form of the ordinary second-order differential 
equations, controlling the behaviour of the gluon profile 
functions in Section [IVJ These equations are highly non- 
linear, and at first sight look fairly untractable. However, 
after introducing a reduction process in Section El we 
can cast these equations in an elegant symmetrical form, 
which allows simple solutions to zeroth order in a s . The 
linear equations for 0(a s )-corrections are also derived. 
In Section [VI] we discuss the (Dirac-like) field equation 
for the quark and the form of the potentials. In Section 
IVIII a self-consistent solution of the coupled quark-gluon 
equations is presented. In Section I Villi we discuss the 
evaluation of expectation values for the energy, spin and 
colour of the quark system. These allow us to define suf- 
ficient boundary conditions to determine the parameters 
of the solutions uniquely. In Section IIXI we discuss the 
necessity to include coupling to GR for the trivial Higgs 
sector. By incorporating vacuum energy as well, one is 
able to make precise estimates of the quark mass and ra- 
dius. In the final section[X]we summarize the significance 
of these results and discuss a possible steps for deriving 
the internal properties and masses of other particles in 
the standard model. 



II. THE QUANTIZATION OF THE EQUATIONS 
OF MOTION IN QUANTUM 
CHROMODYNAMICS (QCD) 

As stated before we start from a basic QCD Lagrangian 
[l8j without electro- weak and Higgs interactions: 

C = i>{i lll D' i -Mty-iF"" .F^, (1) 

where the covariant derivative is defined as follows: 

D^^-ii^A.A". (2) 

Here the A a are the Gell-Mann SU (3)-matrices. The by- 
products run over 8 components: 

8 

A • B = ^ A a B a . (3) 

a=l 

The gluon field tensor is given by the expression: 

F'*" = d"A v - d l, A" + ,g s A^ x A u . (4) 
Here, the SU(3) vector product is defined by: 

8 

(A x B) a = ]T f abc A b B c , (5) 

b,c=l 

where the f a bc are the SU(3) structure constants. The 
classical equations of motion for the gluon fields are given 
by: 

d^(x) = ^g s ij(x)X 7 ^(x) + g s F^(x) x A^x), (6) 
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while the equation of motion for the quark field is: 

(i7^ M + \g s ^\ a A»{x) - m\ ^j(x) = 0. (7) 

In the standard applications of QFT, the fermion and 
boson fields are expanded in the eigenstates of the free 
Hamiltonian, while the interaction Hamiltonian is used 
to calculate the scattering diagrams in the perturbativc 
Feynman series. Its main application lies in scattering 
problems, and the fermion masses are put in by hand. 
In the current application of QFT we start from a more 
basic level where the bare fermions have zero masses (i.e. 
we set M = 0), and we try to solve the full set of field 
equations self-consistently for a particular (quark) state 
vector. The quark is described in its own center-of-mass 
system (the origin of the coordinate system) , and charac- 
terized by localized profile functions representing either 
the gluon fields or internal quark wave function. The ac- 
tual position of the system that describes the effective 
particle is thus irrelevant and only becomes of relevance 
when we enter the quark into the mult i- particle environ- 
ment and consider scattering processes. 

Quantization of the formulation is accomplished by ex- 
panding the fields in terms of creation and annihilation 
operators, the latter satisfying certain (anti-) commuta- 
tion relations. Initially the following expansion of the 
quark fields is used: 

^w = ^w a (x)+^rficw, (8) 

a a 

where b a is the quark annihilation operator, the anti- 
quark creation operator, and the profile functions (j)a{x) 
and 4>a(x) loosely refer to the particle and anti-particle 
wave functions of the quark in its own center-of mass 
system (i.e. the internal wave functions). These wave 
functions are not known beforehand, and even their ex- 
istence is uncertain, since it is not known whether a lo- 
calized system description of the quark state b^O) will 
emerge from the formulation. Since the field equations 
are solved self-consistently, only one space-time coordi- 
nate x enters the considerations, so transitions between 
states are not associated with different times, and the 
overall state vector &|J0) stays the same. Hence, the 
creation and annihilation operators are also time inde- 
pendent and do not evolve over time. The role of these 
creation and annihilation operators in the current for- 
mulation should be distinguished from that in scattering 
problems, where they are often associated with a particu- 
lar fixed functional behavior (a plane wave corresponding 
to a particular positive or negative hyperboloid) . When 
scattering states undergo time evolution through interac- 
tions they change their functional behavior, and creation 
and annihilation operators defined in this way are mixed 
up, allowing even the mixing of creation and annihilation 
operators (Bogoliubov mixing fl9l]). In our formulation 
the creation and annihilation operators refer to the ac- 
tual physical state, represented by frJjO). Isolated parti- 
cles do not change over time (unless they are unstable) 



and their state does not depend on where they are lo- 
cated, so frjj has no space-time dependence. Only, when 
we use an external representation of the particles and 
specify their momentum and/or positions with respect 
to other particles in scattering processes, could we con- 
template time dependent creation and annihilation oper- 
ators to represent the time evolution of the world state 
vector. However, even in that case a time independent 
representation of these operators seems more appropri- 
ate in QFT, as the continuous evolution of the physi- 
cal system between transitions is well represented by the 
continuous profile functions. The evolution of one world 
state to another through transitions is of a stochastic na- 
ture, and should not be represented by changing creation 
and annihilation operators. Rather transition operators 
(e.g. Feynman diagrams) should relate initial and final 
states uniquely defined in terms of fully specified creation 
and annihilation operators, i.e. time-dependent opera- 
tors. Clearly, in our formulation no Bogoliubov mixing 
occurs, and particle and anti-particle components main- 
tain their unique meaning. This is also important for 
the later developments in the theory, when we make an 
important distinction in the algebraic properties of par- 
ticle and anti-particle operators to restore the symmetry 
between these entities which is broken in the Lagrangian 
formulation. In our operator formulation the vacuum 
also maintains its unique identity as the state without 
any particles present. 

We now continue the discussion of Eq.©. The in- 
dex a covers the usual quantum numbers of the standard 
model: spin, color and charge (also called flavor, but we 
only consider one generation of quarks). Since we limit 
ourselves to QCD, we can treat the pair of quark states 
in each generation as a degenerate isospin doublet. The 
family/generation quantum number is supposed to be ex- 
plained by the multiplicity of the solutions and is there- 
fore absent. Both isospin and total angular momentum 
are good quantum numbers in the present formulation, 
so that we can limit the wave function space accordingly 
in the description of the quark state. Since the gluon 
fields carry color and spin, color and spin transitions can 
take place, and an expansion in terms of a complete set of 
color and spin states is required. No expansions in other 
quantum numbers, such as a principal quantum num- 
ber, is appropriate since if a multiple of such states exist 
they should follow as a solution from the given equations, 
rather than appear in the expansion. 

If we use the expansion, Eq. © , in the gluon field equa- 
tion Eq.©, then the gluon field operators also become 
expansions in quark creation and annihilation operators. 
In first instance this expansion looks like: 

a,/3 a,/3 
a,/3 a,/3 

However, the reinsertion of this expression in the Dirac 
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equation, Eq.([7]), leads to higher order terms in the ex- 
pansion of tjj(x), etc., which in turn leads to higher order 
terms in the expansion of the gluon held. It would thus 
seem that we have to revert to perturbative expansions, 
as is done in standard QFT. However, we will show below 
that the operator series can formally be summed to in- 
finity, without increasing the number of profile functions. 

To complete the quantization process we first recall 
the usual anti-commutation rules for the creation and 
annihilation operators: 

= <w, (io) 

and 

|d a ,e^ j = 6 a> p, (11) 

while all other anti-commutators are zero. For free fields, 
when the fermion field is expanded in a complete set of 
plane waves, one can use these anti-commutation rules 
(i.e. their continuous generalizations) to derive the equal- 
time anti-commutator relations. For the current self- 
consistent treatment of the full set of coupled equations, 
such relations are absent. This shows that the anti- 
commutation rules Eqs. (fTU)) and (fTTj) are more funda- 
mental than the equal-time anti-commutator relations. 
Historically one has often assumed the opposite, namely 
that the equal-time anti-commutator relations are fun- 
damental to QFT (e.g. by stressing the analogy with 
the classical Poisson brackets), while Eqs. (flU]) and ([TT]) 
were seen equivalent, but secondary. We have pointed 
out before that it is dangerous to rely on QFT on classi- 
cal analogies, and this is a case in point. 

Since the gluon field operator no longer commutes with 
the quark field operators, it is necessary to specify the 
order of the quark and gluon field operators in accordance 
with the symmetries and physical requirements of the 
theory. The non-commutativity of boson and fermion 
field operators is not new to QFT, for example Bjorkcn 
and Drell [2(| discuss it in the context of QED. However, 
in the current context the consequences are must more 
farreaching. The anti-symmetric nature of the gluon field 
implies that we must replace Eq.Q by: 

•pnv = d^A v -d v A^ + ^- {A" x A u - A" x A M } . (12) 

The quantized gluon field equations can be derived along 
similar lines and read: 

- d^A v + d v (a M A") = |^a 7 > 

+g s {A M x (d^A") - (d^A v ) x A M } 
+|{(9 M A")xA"-A v x(3 fI A") 
-A" x (0"A„) + {d v A^) x A p } 

{A„ x {A v x A p ) - A„ x (A M x A v ) 
-(A"xA")xA )1 + (A f 'xA")xAJ. (13) 



The quantized form of the Dirac equation is given by: 

(i^W - M) if, + ~<? s 7„A a ^ = 0. (14) 

The physical justification for the order of the ip and A£ 
operator will be given in Section I VII 

Before we can derive the full operator expansion of the 
quark and gluon fields we need to refine the algebra of 
creation and annihilation operators. Expressions such as 
ipX^ip display a disturbing asymmetry: the right-hand 
side operator ip contains quark annihilation operators, 
but antiquark creation operators, implying that if such 
an operator is applied to the vacuum, anti-particles can 
be created, but particles can not. However, we know that 
there should be a fundamental symmetry between par- 
ticles and anti-particles and that the decision to call a 
quark a particle or an anti-particle is largely arbitrary. 
So this asymmetric outcome is unacceptable. It does not 
help to invert the role of particles and ant i- particles, as 
then the problem would be that particles can be formed 
from the vacuum and anti-particles can not. It seems 
that the symbolic mathematical language in which we 
formulate QFT is unable to adequately express the nat- 
ural symmetry between particles and anti-particles. The 
solution to this problem is that we can maintain the same 
field theory language, as long as we invert the order of the 
anti-particle operators in the complete string of operators 
between the bra and ket states. This rule can also be seen 
as giving meaning to the physical role of anti-particles as 
being particles moving backwards in time. This elegant 
procedure restores the symmetry between particles and 
anti-particles. The procedure also eliminates the arti- 
fact that QFT gives rise to an enormous vacuum energy, 
which is contradicted by observation. This procedure 
will be indicated by the symbol R, and will be called the 
M-product. Hence, the R-product is a way to introduce 
ordering in expressions that are "timeless" and to en- 
force the property that particles and anti-particles have 
opposite ordering. We should caution, however, that this 
R-product only applies to quantum operators with iden- 
tical space-time variable x (as is the case in our current 
application). Operators referring to different variables x 
and x' are not subject to this product, as their time-order 
can already be specified by the different times t and t' . To 
ensure the correct application of the R-product one could 
append a variable label to the operators, so that oper- 
ators referring to different space-time variables are not 
confused with those having common variables. Since the 
product does not apply to operators at different points, 
it has no effect on most scattering series in QFT, which 
may explain why it has not been discovered before in 
the application of QFT. The product is however already 
used implicitly in standard applications of QFT, when 
the normal product is inserted to remove unphysical vac- 
uum terms. We actually discovered this rule numerically 
before its theoretical origin was uncovered, as the coupled 
equations required unexplained minus signs, which could 
only be justified after the R-product was introduced. The 
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detailed explanation of this product is given in Appendix 

HI 

Using this ]R-product one can derive a formal solution 
of the quantized equations of motion (Appendix [Bj . In- 
stead of the expansions, Eqs. flSJ and ([9]), we get the 
exact expressions: 

if>(x) = Soo ]T {b Q ^ a (x) + 4c (*)} , (15) 



and 

+ 2 ilSoodjA^ + £ d-Sooft^A^, (16) 

ct,P a,P 

where is the result of taking the operator 

2-N-N n-N-N 



S„ = (l-iV-iV> 



2 n 
n — N — N 
n 



(17) 



to the limit n — > oo. The operators -/V and JV are given 
by: 



N = Y,b%, 



and 



A 



(18) 



(19) 



The operator Eoo can be recognized as vacuum projection 
operator, as it only gives non-zero (unity) results if N — 
N = 0. As far as we know this operator has never before 
been encountered in QFT. It will play an important role 
in the solution of the self-consistent field equations. Since 
this operator is surrounded by creation and annihilation 
operators, it effectively acts like a one-body, rather than 
a vacuum, projection operator. 

The detailed derivation of these relationships is given 
in Appendix [B] Since the space of internal states for the 
quark is finite (see Section IIIip , this infinite expansion 
of Sqo is terminated in practice after a finite number of 
terms. Clearly, this is not true if the state is characterized 
by a continuous degree of freedom, such as in the case 
of scattering, when we have to define N as an integral 
over all momentum states. An important property of this 
expansion is that the original number of profile functions 
in the approximate expressions, Eqs. (JSJ and ©, does 
not increase in the transition to the exact expressions 
(fT5j) and (fT6|) . This property is largely responsible for 
the practical feasibility and the exactness of the current 
theory, and ultimately for the power of QFT. 

We now want to illustrate the effective one-body nature 
of the projection operator. If one operates with ij>(x) on 
a one-body state one gets: 



however, for a many-body state one gets a zero result: 

r/>(x) 1*4- ••&+„! 0)=0, n>l, (21) 

and similarly for anti-particle states. This shows that 
the only fermionic solutions of the coupled equations in 
the current formulation are of one-particle nature (the 
term particle is used here as a generic term for particle 
or anti-particle). This clearly demonstrates that we can- 
not apply this formulation to a composite problem like 
a proton system, as this involves at least three quarks. 
The only way to generalize this approach to composite 
systems is to collapse the three quarks into a single color- 
singlet point-particle. However, this would eliminate the 
QCD interaction and thus would not yield a proper de- 
scription of nucleons. For leptons this might be a possible 
model, however, this would require the introduction of 
the electro-weak interaction at this basic level, which we 
have not attempted yet. Obviously, for protons and other 
hadronic systems the standard lattice gauge calculations 
remain the best treatment [2l| . Now that we have solved 
the quantum operator problem formally, we can proceed 
towards the determination of the profile functions. We 
will see that QFT has more surprises in store, and that 
the simplest solution for the single quark system is of 
great beauty, and realizes the localization of the bound 
system in an amazing way, never seen before in physical 
theories. 



III. EXPRESSIONS FOR THE PROFILE 
FUNCTIONS 

The next step after having removed the (creation- and 
annihilation-) operator structure in the field equations is 
to remove the color and spin dependence, by introducing 
scalar wave function components and scalar gluon profile 
functions. We start with the paramctrization of the wave 
function of the lS-quark state: 



b a (r,t) = exp(-iEt) 



1 



ier • rg(r) 



a • 



where \a is the spin wave function, rj a is the isospin 
wave function, and £ a is the colour wave function. Here 
we assumed that the radial wave functions / and g are in- 
dependent of spin, isospin and colour. This charge inde- 
pendence is no longer valid if electromagnetic forces are 
introduced, and u- and rf-quarks become distinct. The 
corresponding anti-quark states is given by: 



C(r,t) = exp(LEi)- 



1 



ier • rg(r) 

m 



XaVaZa- (23) 



^(x)\bi\0) = (P a (x)\0), 



(20) 



Notice the absence of phase factors. Usually, one defines 
the spin function for the anti-particle by expressions such 
as (— ic2Xa)- Also, the colour and isospin wave functions 
for the anti-particle belong to adjoint representations of 
the relevant SU(3) and SU(2) representations. However, 
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the particle and anti-particle creation and annihilation 
operators automatically take care of these phase factors. 
Hence, Eq. (|23|) does not correspond to the normal defi- 
nition of the anti-particle wave function, but exploits the 
simplification enabled by the operator formalism to ex- 
press the anti-particle wave functions in the same wave 
function basis as the particle states. The time depen- 
dence of these wave functions reflects the fact that we 
describe stationary states. Although, the quarks belong- 
ing to the higher generations will decay into the lighter 
quarks via several decay mechanisms, we expect that 
these higher generations are also characterized by a real 
positive eigenvalue E, and that the decay width is deter- 
mined in the scattering problem and not in the current 
bound-state formalism. 

The wave functions / and g are normalized according 
to: 



dr[f(r) 2 +g(r) 2 ] = l. 



(24) 



IV. DERIVATION OF THE SCALAR GLUON 
FIELD EQUATIONS 

After eliminating the common quantum operators one 
is left with equations in terms of the remaining SU{n) 
operators and functions. The SU (3) operators can easily 
be factored out using the relationship: 



E 

b,c=l 



fabc^b^c — 3iA a . 



(31) 



The spin operators satisfy numerous spin identities, 
which can be used to remove their explicit occurrence 
as well. We are left with a set of scalar equations of 
motion for the gluon profile functions. To simplify our 
notation we introduce some short-hand notations for the 
source functions. We define: 



S 



as P+9 2 



(32) 



Using the wave functions Eq. (|22"j) and Eq. ([23|) . we can 
now evaluate the source terms in Eq. (fl"3|) . and thereby 
determine the form of the different gluon field compo- 
nents in Eq. ([TB|) . We express the results in terms of five 
profile functions: 



aO, pp 



(x) 



Fo(r) 



rspin /\ \ 
%3 V A a)af: 



(25) 



\PP 



(x) = -A 



0,aa / 



F(r) 



(r x (r) afj (\ a ) a0 , (26) 



.F (r) fA 



= • °U (Aa) Q0 exp(2iM), (27) 



^{x) = -i^{i'<rU{K) afi exp(-2LBt), (28) 



a pa I \ j F i{r) i F 2 (r) 
L ys 1 ys < 

x ( X a) a p exp(2i£t), 



K P aff( X ) = { 



Fi(r) 



-Ca/3 



F 2 (r) 

9sT g s r 
x (A Q ) Q/3 exp(-2LEt). 



f (f • <T) a (: 



(29) 



(30) 



The dummy isospin Kronecker 5-function has been sup- 
pressed throughout. This simplification becomes clearly 
invalid as soon as we introduce isospin breaking QED 
forces. In the next section we will derive equations for 
the five profile functions F (r), F(r), F (r), i*\(r) and 
F 2 (r). The Dirac equation yields equations for the two 
functions / and g, and are discussed in Section [VII 



c a s 2fg 
2 r 



S 2 



a s f 2 -g 2 



(33) 



(34) 



where a s is the strong coupling constant a s — g 2 /Air. 
These source terms are related by the identity: 



q2 q2 q2 

°0 ~ °2 — °1 • 



(35) 



We also introduce the auxiliary profile function K , which 
is going to play a central role in the dynamics: 



K = F 1 +F 2 . 
The differential equation for F reads: 



(36) 



F" =S - 4^1(1 - 3F) - \(2Er + 3F Q ) 



x (2F{ 





2 , R 2 



K(rF' - F ) 



3 ^ 



F (rK' + K). (37) 



As expected, quadratic and cubic terms play an impor- 
tant role in these equations. These terms express the 
non-linear nature of QCD, and the self-coupling of glu- 
ons. The solution of the linearized equation would lead 
to a Coulomb like potential. 

For the magnetic particle-particle component F we ob- 
tain: 



F" = S 1 -^F F 1 (2Er + 3F ) 
|F(1 - §F)(1 - 3F) + A(F 2 - F*)(3F - 1) 



%K(rF[ - Fi) + \f x (tK> + K) - ^FF X K. (38) 
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An earlier version of Eqs. (j37|) and (|38|) . without particle- 
anti-particle contributions, was presented in Ref . [l6| . 
These equations did not lead to bag-like bound states. 
However they already displayed the very interesting non- 
linear behaviour, suggesting soliton-like solutions. It has 
since become clear that the particle-anti-particle coupling 
is essential for the appearance of singular confining po- 
tentials. 

The profile function Fo, corresponding to the gluon 

ip , 



field operator A^^gix), satisfies the equation: 



and 



G = 3F 4 2 - ZFf - 



1 



We then can write Eqs. (|37H4T)|) as follows: 



F a — S 



6 



7F1H 



; F 3 K 2 - -KF> 



^-F (rK'-K), 



(45) 



(46) 



F o =Si+ (1 - 3FY - -K 



~^K(rF^ - Fo) - ^(2ER + 3F Q )(rK' + K) 

+ ^F 1 (l-3F)(2Er + 3F ). (39) 

Finally, the two profile functions Fx and F 2 appearing in 
A^ a a/3 (x), satisfy the equations: 



F% =S 1 - ^F 4 H - ^F K 2 - - r KF> 
~Fi{TK'-K), 



F'i 



Si tjFoH H — 2GF4 

-*KF[- 1 F x {rK' -K), 



(47) 



(48) 



F» = S + ±(l-3F)(rK>-K)-^ 
^(2Er + 3F ) 2 - ^F (2ER + 3F )(1 - 3F) 



1 



- 3F) 2 



F'K 
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F X K 2 



(40) 



and 



F' 
1 2 



3^2 



S2 — So 



Fo 



{2Er + 3F ) 



-F{' - Mi + 5£i + *pa K _ {2Er + 3Fo 



F' 



F 



Fa 



(2Er + 3F )(l-3F) 



rpt Fo 



3F 
3 



F(F^ + 3F[) 



■ —F'F 2 — ~^ 2 ~FF 2 — ~^FF\ 
^-F 2 (2K-F 1 ) + ^F 1 F 2 . 



(41) 



These 5 equations for the 5 profile functions do not look 
very tractable, however, after introducing a number of 
new auxiliary functions they take on a much more elegant 
form. We introduce the alternative profile functions: 



^3=^0 + ^ Er 



(42) 



and 



F{> 



So - ^F 3 H + ^GFi - ^FiK 2 
-- r KF' A -^F A {rK'-K). 



(49) 



This has not just simplified the equations enormously, 
but also has brought out a strong symmetry between F3 
and Fo on the one hand, and F\ and F4. on the other. 
The equation for F 2 , Eq.lflT]). can be added to Eq. (|30|) . 
to yield a very elegant algebraic expression for K: 



K = r 



■S 2 - 3F 3 F^ + 3F^F - 6F4F1 + 6F{F 4 



9FI 



9F 2 



18F 2 



18F| 



(50) 



We can derive another important relationship by differ- 
entiating the numerator in Eq. (|5"0"|) . and inserting all the 
expressions for the second derivatives. We obtain: 



1 



S (F + 2F 4 ) = 5i(F 3 + 2Fi) - -(rS 2 y 



(51) 



If we combine this equation with the Dirac equation, we 
can derive a simple relationship for the source functions, 
which is related to current conservation. In the next sec- 
tion we continue our treatment of the gluon profile func- 
tions exploiting the symmetry apparent in the equations 
(|46H49p . These steps are necessary to progress towards 
a physical solution, however, despite their elegance these 
steps are of necessity rather technical. Hence, those who 
want to immediately progress to the treatment of the 
quark field equations can skip to Section IVII 



Fa 



(43) 



Furthermore, it is convenient to introduce the following 
expressions: 



H — 3F3F1 — 3F4F0 



(44) 



V. REDUCTION OF THE SCALAR GLUON 
FIELD EQUATIONS 

We can exploit the symmetry of the equations (|46I 
to partly decouple the equations. To this end we 
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introduce further auxiliary functions: 

X = F 3 — F ; Y = F 3 + F ; 
U = F i -F 1 ; V = F 4 + F 1 ; Z = K/r. (52) 

By adding and subtracting the equations (|46H49[) we find: 
X" = S- + ^U(VX-UY)-9Z 2 X+6ZX'+3XZ', (53) 



Y" = S+--2 V(VX-UY)-9Z 2 Y -6ZY' -3YZ' , (54) 



U" = -S- + ^X(VX UY) + ^(9UV 1) 



-9Z 2 U + 6ZU' + WZ', (55) 



V" = ~ £2* (™ - UY) + J (9C/T> - l) , (64) 



where 



S+ — S+E; S— = S—E 



(65) 



We can also express Eq. ([501) in terms of these new func- 
tions, and obtain: 



1 



XY 1 - YX' + U'V - UV' 



1 



:rSo. 



(66) 



We can easily solve Eqs. (l6H6"4f to 0(a°), by ignoring the 
source terms. By imposing the parity requirements of the 
potentials (see Section IVI[) we are lead to the following 
solutions: 



X = Y = j3Er; U = V 



0, ±- 



(67) 



V" = ^-Ay{VX-UY) + ^(9UV-1) 



where 



2r 2 ' ' 

-9Z 2 V - 6ZV' - WZ', (56) 



S+ — So + Si] S— — So — S\. (57) 



We have partly decoupled the equations and can intro- 
duce further simplifications by introducing the expres- 
sion: 

/ r 

Fin = rxp [ :-! / </r'-^— I = exp I 3 / <l r' Z:.r' ■ 

•(58) 

Because of the symmetries of the equations one can also 
define E with a minus sign. We can invert this relation 
by expressing K in terms of E: 



K{r) 



r E'(r) 
3 E(r) ' 



(59) 



Again this result is insensitive to the sign of E. We now 
write: 

X = EX; Y = E- 1 Y; U = EU; V = E~ 1 V.(60) 

The new functions X, Y, U and V satisfy even sim- 
pler coupled equations, as the trailing Z dependent terms 
have vanished: 



X" = S- + (\'X ry 



Y" = S 4 



VX - UY 



U" = -S- + ^X (VX -UY)+^ (WV - l) , (63) 



(61) 



(62) 



The solutions for a — ^ and a = — i appear to be phys- 
ically equivalent, the former corresponding to the nega- 
tive sign of F. Hence, we will only consider the set with 
a = —1/3. We will see in Section lyilll that the a = 
solutions lead to infinite energy integrals, and must be 
excluded. For a = — 1/3 the |/3| = 2/3 solutions are of 
special interest. For j3 — 2/3 we obtain a trivial solu- 
tion where all potentials are zero (Fo = F = Fq = Fi = 
0; K = 0; E = 1). For /3 = -2/3, K - and thus E - is 
singular at a finite radius tq (the bag solution) . As a con- 
sequence, all original profile functions F3, Fo, F4 and Fi 
are also singular at ro- However, the reduced functions 
X, Y , U and V remain finite. Hence, our reduction of the 
equations in this section does not just lead to a simplifi- 
cation of the equations, but it also leads to a systematic 
elimination of the singularities from the equations. 

The solutions summarized in Eq. (l67[) are valid in the 
absence of source terms. The existence of solutions for 
a s J. is an important consequence of the non-linear 
nature of the field equations in the current bound-state 
approach. It is an illustration of the power of non-linear 
Gauge theories in the construction of emergent complex 
structures. In the standard perturbative approach to 
QFT, a theory with a s 4- would make no sense, and 
there certainly would be no binding. 

Although terms higher order in a s do not play an 
instrumental role in the binding mechanism, they will 
modify the detailed outcome of the theory. In the fol- 
lowing we show how these terms can be treated per- 
turbatively, expressing the corrections in terms of lin- 
ear equations for the remaining profile functions. Com- 
bined with the calculation of observables in Section [Villi 
these linear equations provide a powerful tool in analyz- 
ing the physical solutions in more detail. The combina- 
tion of exact 0(a°)-solutions, with perturbative Over- 
corrections might look somewhat heuristic, however, non- 
linear equations require novel methods, and this is a very 
effective approach The natural place to start the pertur- 
bative scheme is the set of regularized equations Eqs. dGlt - 
I64p . which are defined in terms of non-singular hatted 
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profile functions. However, the functions X, Y, U and 
V have non-zero base values. For a perturbative scheme 
it is preferable to start from zero functions. This can be 
done by reverting back to functions like Fq, etc., but now 
defined in terms of the hatted functions. Unfortunately, 
this requires a new set of auxiliary functions. However, 
the symmetry will not escape the reader, and the scheme 
remains very elegant and natural. The new function are: 



These equations can be solved formally in terms of the 
source functions. For the two uncoupled equations we 
obtain: 



and 



F Q = X r - r / dr'So - / dr' r' S 



(79) 



F 3 = {X + Y)/2; F Q 

F 4 = (U + V)/2; F 1 

F = F 3 -(3Er; F 

H = 3F3F1 — 3i^4-Fo! G 



= {Y-X)/2; 
= {V-U)/2; 
= i*4 - a; 

= 3F| - 3ff - i (68) 



For P = 2/3 and a = -1/3 one gets K = and E = 1, 
and the old functions are recovered (i.e. Fq — > Fq. Ex- 
pressing the differential equations into the new auxil- 
iary functions, we are led back to differential equations 
of the familiar form (|46M9[) . but now without the K- 
contributions: 



pn _ c 6 
^3 - *0 - ^ 



Fq = S\ yF 4 H, 



32. 3 „ 
F A = S\ tjFqH H — ^Fi G, 



F" — So 



-Ft G, 



(69) 



(70) 



(71) 



(72) 



where the reduced source functions are defined by: 

S = (S+ + S-)/2; Si = (S+-S-)/2. (73) 

In terms of these new functions the 0(a!?)-solutions 
Eg. (1571) are given by: 

F 3 = /3Er; F 4 = a; F = A = F Q = F = 0. (74) 

where a — 0,±i as before. We now expand the differ- 
ential equations to first order, setting a = —1/3. The 
0(a s ) linearized equations now read: 



F — So, 



6/3E 



2 ~ 



F = Si H F\ H — tjFq 



F{' = S - 9p 2 E 2 F 1 - ^Fq 



(75) 
(76) 

(77) 

(78) 



fir 



dr 



1 

3^ 



dr 1 r' 2 Si, 



(80) 



where we assumed that the physical solutions have a fi- 
nite spatial extent specified by tq. The parameter Ao 
could be combined with (3 into a single parameter. How- 
ever, we prefer to specify the 0(a° )-solution by a fixed 
value of /3, and treat the 0(a s )-corrections via Ao- The 
other two equations can be decoupled. First we multiply 
Eg. ([76]) with r, and then differentiate the result twice. 
After eliminating the function Pi, we can express the 
result as follows: 



P"{r) + -P'(r) + 9f3 2 E 2 P(r) = r" 2 Q(r), 
r 

where P(r) is defined by: 
P(r) 



Fo 
r 



(81) 



(82) 



and Q(r) is defined in terms of source functions by its 
derivative: 



Y(r) = (rSi)" + 6/3ES + 9f3 2 E 2 rS 1 



(83) 



One can re-express this in the original source functions, 
and show that: 



Q(r) = (rSi)' + 3/3ErS 2 + u>, 



(84) 



where oj is the third integration constant. The solution of 
Eq. (f8"T)) can be expressed in terms of the spherical Bcsscl 
functions ji and ni: 



P(r) 



>Ji{ar) 



v- dr'r'ni{ar') Q(r') 



a2 ni(ar) / drVji(ar/) Q(r/)j (85) 



where a = 3f3E and v is the fourth - and final - free 
constant of the 0(a s )-solutions. For the special case 
\fi\ = f 1 Q = w to 0(a s ). For small r a non-zero ui 
generates a singular component ln(r) in P(r) and a com- 
ponent r 2 ln(r) — r 2 in Fo . The other function Fi can be 
deduced directly from Eo. fTB")) : 



2aF[ = -(rSi)' ~ a 2 r 2 P{r) + Q(r) = 
= -a 2 r 2 P(r) + arS 2 (r) + w. 



(86) 
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We can get one expression without derivatives from 
Eq.®: 



- + aF 1 = --(r5 1 ) + - 



P'ir) + -P(r) 
r 



(87) 



At this stage we have four free parameters: Ao, ^, v 
and w, with dimensions E, E 2 , 1 and -E 2 , respectively. 
These parameters should all be of 0(a s ). The imposition 
of such a constraint on " free" parameters seems contra- 
dictory, but the non-linear nature of the full equations 
limits such freedom in the full solution. We can enforce 
the 0(a s ) constraint by absorbing the free parameters 
into the left- most integrals in Eqs.([79]l. ([80]) and <j85j), by 
allowing the upper integration boundary, which is now 
fixed at ro, to vary between and tq. In this way the 
corresponding parameters are always of 0(a s ), and also 
have the right order of magnitude. The need for such 
a prescription is a consequence of the hybrid nature of 
the solution scheme: one tries to built a perturbative so- 
lution on a non-linear O(ttg) basic solution, which itself 
contains a free parameter f3. 

This completes the discussion of the equations for 
the gluon profile functions. Although we seem to have 
strayed a long way from the original gluon field tensor 
defined in Eq.Q, expressing this tensor in terms of the 
reduced profile functions is straightforward. In the next 
section we discuss the Dirac field equation, which will 
give further constraints on the profile and source func- 
tions. We stress again that all profile functions and quark 
wave function components have to be determined self- 
consistently. Hence, although we made great strides in 
deriving simple equations for the gluon profile functions 
in this section, we need the Dirac equation to complete 
the solution. 



VI. THE QUANTIZED DIRAC EQUATION FOR 
ISOLATED QUARKS 

We have to quantize the classical Dirac equation, 
Eq.©. Since both ip and A% are operators expressed 
in quark creation and annihilation operators, quantiza- 
tion amounts to deciding in which order these operators 
have to appear in the quantized equation. As stated in 
Section [Til the correct choice is ipA% , leading to equation 
(UJ|. This choice guarantees that the quark interacts 
with the fields it generates itself. The opposite choice, 
A%ip, does not yield a bound-state solution. For exam- 
ple, it gives rise to the operator b^bpb^, which vanishes 
when operating on a one-particle state bg |0). Since two- 
and higher-body terms are eliminated by the "projec- 
tion" operator Sqo, none of the other terms survive ei- 
ther. On the other hand, ipA£, leads to the operator 
b^b^bp, which yields a finite result, when operating on 
a one-particle state b\ |0). Hence, the only way we can 
get binding potentials and bound normalized states is to 
choose the indicated order. Again we have to extract 



sets of terms from the equations with similar operator 
form. We then can eliminate these operators and are left 
with algebraic equations. In the reduction of the opera- 
tor equations we have to take into account the M-product 
when manipulating the anti-particle terms. Ignoring this 
aspect would lead to erroneous signs in the reduced equa- 
tions and inconsistencies. 

We can express the reduced Dirac equations as a set 
of scalar equations for the large (/) and small (g) wave 
function components. We find: 



(E - V - M - V.) f(r) =-i^±-g>(r) + V T g(r), 



and 



{E-V + M + V s ) g(r) 



f(r) 



+ f'(r) + V T f(r). (89) 



The bare quark mass M is zero in our application, how- 
ever, we leave it here as it could be non-zero when a non- 
zero Higgs field is introduced. Our aim is to derive the 
masses of dressed quarks, so it would be inappropriate to 
introduce phenomenological masses at this point. 

The formulation gives rise to three types of potentials 
(if we introduce isospin breaking effects we would get a 
fourth type of potential, as is allowed by the symmetries 
of the Dirac equation. This fourth class is also of tensor 
type). The vector potential is given by: 



V(r)=--(F +2F 1 ). 
r 

while the tensor potential is given by 

C 



V T {r) 



(f 



2F 



(90) 



(91) 



Finally, the scalar potential is exclusively given in terms 
of Kir) : 



V s {r) = —K{r). 
r 



(92) 



In these expressions C is an SU(3) constant with the 
value 8/3. Each of the potentials receives contributions 
from the particle-anti-particle gluon profile functions. 
The latter terms are essential for the emerging elegant 
self-consistent solutions of the equations and thus indis- 
pensable for the existence of dressed quarks that are fi- 
nite in extent. Although the particle- anti-particle gluon 
amplitudes fluctuate rapidly with time (they are propor- 
tional to exp2i£'i or exp— 2iEi), this time dependence is 
cancelled out in the scalar Dirac equations. This is one of 
the many elegant features of the current self-consistent 
theory. The required behaviour of the potentials near 
r = fixes the parametrization in Eg. (1671) . In lowest 
order 0(a°) we demand that V and V s (which are even 
in r) are constant near r = 0, while the tensor potential 
(which is odd in r) behaves linearly in r for small r. This 
is consistent with the small r-behaviour of the odd wave 
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function f(r) ~ T and the even wave function g(r) ~ r 2 . 
We will see later that in 0(a s ) the lower component g(r) 
also has small r components r 2 ln(|r|). 

By multiplying Eq.([55|l with g(r), and Eq.® by /(r), 
and adding them together, we can derive the following 
relationship: 

(rS 2 )' = 2(E - V)Sir - 2rV T S a + 2S a . (93) 

Combining this with Eq.(|51j) we obtain: 

S V T + S 1 V = 0, (94) 

and 



2S + 2ErS 1 = (rS 2 )'. 



(95) 



Eq. (|95[) also follows immediately from current conserva- 
tion: 



d» = 0, 



(96) 



if we take the particle component on the right, and the 
anti-particle component on the left. Notice, that Eq. (j9"6")l 
is trivially satisfied for the particle-particle components. 
Hence, our derivations so far have been completely con- 
sistent with current conservation. We can derive another 
equation from Eqs. (|88H8"9l in terms of source functions: 



(E - V)rS 2 = rS (M + V s ) - ^(rS^ '. (97) 

We now derive the basic non-linear self-consistent solu- 
tion of the QCD field equations by combining the condi- 
tions on the gluon and quark profile functions. 



VII. THE BASIC SELF-CONSISTENT 
SOLUTIONS FOR THE DRESSED QUARK 

We anticipate that if a bound solution exists for the 
dressed quark system, that it is bound by strong, possibly 
singular forces to ensure that its size is small enough to 
be consistent with the known properties of quarks. The 
profile function K may well play a prominent role in this 
regard, as it is given by a ratio (see Eq. (|50"|) ) and thus 
would become singular if the denominator turns zero. 
The function K stands alone amongst the other profile 
functions in being given by an algebraic form. This form 
is a unique consequence of the non-linear nature of the 
theory and the inclusion of both particle and anti-particle 

configurations. The reduced profile functions Fq, F, Fq 
and Fi satisfy regular second-order differential equations 
(|75H78j) . and are unlikely to have singular or bag-like char- 
acteristics (we will show in the following that the reduced 
source functions in these equations are also devoid of sin- 
gularities). We should note that singularities in K or E 
would not just be reflected in V s , but also in the other 
potentials, because these potentials are implicitly defined 
in terms of E and 



To analyze the behaviour of K or E we cannot di- 
rectly use Eq. (f50| . as we have no simple equations for 
the unreduced profile functions. In addition, the unre- 
duced profile functions may well be singular. Rather, we 
start by casting Eq.(f5Tj) in reduced form: 



-S (F +2F i ) + S 1 (F 3 + 2F 1 ) 



(98) 



The fact that the new equation has the same structure 
as the unreduced equation illustrates the elegance of the 
reduction process. The advantage of the reduced form 
is that it is expressed in the regularized profile functions 
with their known regular 0(a°) and 0(a s )-solutions. By 
re-expressing Eq. (|98j) in terms of the original physical 
source functions, which must be normalizable and display 
regular properties, we can derive explicit expressions for 
E and E' 1 ; 



(rS 2 y±yJ[(rS 2 y] + S 2 (A 2 -5 2 ) 
E = - — — — , (99) 



and 



E- 1 
where 



S+ (A — B) 
<rS 2 )' ± ^[{rS 2 )'} 2 + Sl{Ai-B*) 



S- {A + B) 

A = 3F 3 + 6Fi -> 3/3Er, 
B = 3F + 6F 4 -> 6a -> -2. 



(100) 



(101) 



Here the right-most expressions are the 0(a®) solutions. 
We can use Eq. ([93|) to rewrite the derivative term. 

For p = 2/3 and the top sign in Eas. ([99lH00]l we obtain 
E = 1 and K — 0, corresponding to the trivial solution 
with zero potentials. The mirror solution, with /3 = —2/3 
and the lower (negative) sign in Eqs. (pM UOOp . yields the 
non-trivial bag solution. For both these solutions the po- 
tential singularities of E or E~ l , which would be due to 
a zero of 3(3Er ± 6a, cancel out. Notice also that the 
reduced source functions remain finite, even when E or 
E^ 1 become singular when S+ or S- passes through zero. 
Hence, the reduced gluon Eqs. (|75H75|) are devoid of sin- 
gularities, and form a good starting point for calculating 
higher order effects. 

Taking j3 = —2/3 we have: 



(/ 



(f + gy 



(102) 



Here E is singular for / = —g, which serves to define 
the radius of the system. This condition coincides with 
the usual demand in the MIT bag model [22J at the bag 
surface. Using Eq. (|59l) we obtain: 



K 



4 f'g-fg' 
3 r p-g 2 ■ 



(103) 
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This result can also be obtained by explicit calculation 
from Eg. ([5U| . Hence, the scalar potential equals: 



V s 



4C f'g - fg' 
3 f 2 -g 2 ' 



(104) 



where C = 8/3 as before. Using Eq. (ll02[) we can also 
derive the other 0(a°)-potentials: 



V = 



AC f 



3r (p 



and 



8C fg 



3r (f 2 



[Er(f 2 +g 2 ) + 2fg] 7 (105) 



[Er(f 2 +g 2 )+2fg] . (106) 



If r = 7"o is the point where / = —g, then the quark 
system is confined to r < r$, as the potentials become 
infinite at Tq. Hence, if this solution exists, then it is a 
bag solution with confining potentials Eqs. (jl04fll06|) . 

Can we construct a solution of the Dirac equations for 
the potentials Eqs. (|lU4H106[) ? Clearly, this is a highly 
self-consistent problem, as the potentials are expressed 
in the wave functions. We consider the Dirac equations 
Eqs. (|88H89|) in more detail. First, we combine the vector 
and tensor terms. Because of the identity, Eq. ([M|) . we 
can rewrite these equations as: 



and 



(E-M-V a e ")f(r) 



(E + M + V s eff ) g{r) = 



g( r ) 



g'(r), 



+ f'(r). 



(107) 



(108) 



with the potentials Eqs. (1104H106|) . The solution is given 
by: 



f(r) = Nsm(Er) 



and 



g(r) = N 



cos(Er) 



sin Er 



Er 



The wave functions are defined in the range: 

Q < Er < Er = x = 2.04278694 
leading to a normalization constant of: 



N~ 2 = 



Xo 



sin x \ ! 

x 



(111) 



(112) 



(113) 



(114) 



These solutions allow us to calculate the scalar, vector 
and tensor potentials. In Fig. [I] we display these poten- 
tials. The origin of the length and energy units used in 
this figure will be explained in Section IIXI Clearly, the 
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where 



V. 



eff 



V S + V^ = V S -V T ^. (109) 



V^f is a convenient mathematical vehicle for construct- 
ing a solution. Notice that Eqs. (|107H109p are completely 
general and do not depend on the specific 0{oP s ) solution 
considered in Eq. (|101j) . 

Going back to the specific solution, we insert the po- 
tentials Eqs. (|104H106)) into Eq. (fT09l . We then get: 



4C rjf'g 
3r 



fg')-Er(f 2 +g 2 )-2fg 



f 2 



(110) 



We can now eliminate /' and g' in this expression by 
means of the effective equations Eqs. (|107p and (|108l) . For 
M = one then finds V s eff = (AC/ 3) x V s eff , which only 
can be satisfied for — as C = 8/3. This yields the 
free / and g solutions, except that they are now confined 
to the volume up to the singularity in the potentials. Re- 
inserting these solutions back in Vf™, one finds that this 
effective potential is indeed zero. Hence, we have found a 
self-consistent solution of the Dirac equations for M = 0, 



r[fm] 



FIG. 1: Confining potentials for a u or d quark 

vector and tensor potentials are more singular than the 
scalar one. In Fig. [5] we show the resulting wave func- 
tions, Eqs. (|HHI114l) . normalized over the interval [0, ro]. 
These wave functions (in different units) also appear in 
the phenomenological MIT quark bag model j22|, where 
they describe the wave functions of the three quarks in 
a nucleon. Hence, the current theory shows how such 
a bag model could possibly originate from fundamental 
QFT considerations, although the restriction of our the- 
ory to single-particle states forbids a direct application 
of our theory. Obviously, the context is also entirely dif- 
ferent as the size of an MIT bag is about 1 fm, while our 
single quark bag has Planck length dimensions. 

We derived the solutions for a zero bare quark mass. 
If we introduce a non-zero Higgs field the mass M in 
the Dirac equation could be non-zero. In this case we 
find that is no longer zero but becomes = 

— (32/23)Af. Hence the addition of a bare quark mass 
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r[fm] 



FIG. 2: Quark wave functions 



we will also derive the formula for the total (internal) en- 
ergy of a quark system. The zeroth order result will turn 
out to be negative. This presents a serious problem if we 
want to identify this energy with the effective mass of the 
dressed quark. However, there are various reasons why 
the model cannot be complete. In the absence of the 
Higgs field there is no scale parameter in QCD, hence 
the scale of the quark system is undetermined and it is 
a mystery how the small quark mass of the light quarks 
can be reconciled with the effective pointlikc nature of 
the quarks. In Section IIXI we suggest a novel way to 
complete the model of single quarks leading to a solution 
of these interconnected problems. 



M leads to a solution with an apparent mass of M — 
(32/23)M = -(9/23)M. This is an illustration of the im- 
portance of self-consistency in these bound-state QFT so- 
lutions, and shows that any perturbative considerations 
may fail, since they even would suggest the wrong sign 
of the mass. It shows that the calculation of quark solu- 
tions for non-zero Higgs fields is feasible, although they 
have important further consequences, such as additional 
terms in the total energy expression. 

The solutions obtained in this section also have to sat- 
isfy additional consistency constraints. If we describe an 
isolated quark by the state b£ o) , then its isospin, color 
and spin are specified by the index a. However, these 
properties can also be defined as expectation values and 
expressed as volume integrals over the quantum fields. 
Naturally, the two possibilities should yield the same 
answer. For isospin this condition is trivially satisfied, 
however, for spin and color things are more complicated. 
These conditions will be discussed in Section FVlIIl Our 
0(a°) solution for (3 = —2/3 satisfies the resulting con- 
ditions exactly. Hence,we have found an exact solution 
of the equations of motion for a s — 0. From the analysis 
of scattering problems we know that QCD is an asymp- 
totically free theory. Hence, it may be quite reasonable 
to set a s — at the Planck scale. Since this is the scale 
that seems to be relevant for the light quarks (see Sec- 
tion |IX| the a s — 0-case might a complete description of 
dressed quarks in the QCD context. However, in order to 
describe the mass splitting between the quarks of differ- 
ent charge we need to introduce QED. This theory is not 
asymptotically free, so that we cannot ignore this cou- 
pling and have to apply a perturbative approach in terms 
of non-zero coupling constants at the Planck scale, any- 
way. If we assume that a s ^ 0, then we have to consider 
higher order corrections, as well and impose the bound- 
ary conditions discussed in I Villi The color condition is 
the most stringent one and leads to three separate condi- 
tions in 0(a s ). Hence, it effectively fixes the integration 
constants v and u, while putting indirect conditions on 
Ao, which in turn leads to 0(a s ) corrections to the value 
of the radius Xq. The value of fi is essentially fixed by 
the spin condition. Values for v and w for the bag so- 
lution are given in the next sectior fVIIIl In that section 



VIII. OBSERVABLES AND CONSERVED 
QUANTITIES 

The symmetries of QFT determine the conserved quan- 
tities and observables. The invariance under phase trans- 
formations of the quark field operator leads to current 
conservation (cf. Eq. ([M| V This relation implies the con- 
servation of baryon number B: 



B = ~ I d 3 x B n ). 



(115) 



where the state vector is either a quark: 

\B a ) = bt |0) , (116) 

or an anti-quark: 

\B a ) = 4 |0) . (117) 



The normalization in Eq. (|115l) is chosen in such a way 
that the quark has baryon number 1/3, and the anti- 
quark — 1/3. Notice that we need the R-product to obtain 
the latter result. 

The invariance under translations leads to the canoni- 
cal energy-momentum tensor: 



_dC_ chp_ _d^dA% _ 



(118) 



By adding a complete differential one can re-express this 
as the Gauge invariant expression ((23[): 



&» = -CgT + fofD^ - F^F^ p , (119) 



where: 



d v T^ v = cW 1 ' = 0, 



(120) 



so that the spatial integrals of 0^° are conserved quan- 
tities, corresponding to total energy for \i = and total 
momentum for fi = i. It is easy to quantize Eq. (|119p by 
taking A£ outside the quark current, symmetrizing the 
gluon term, and introducing the R-product. 
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The spin of the quark is related to the tensor: 

M^ p = M^ p + d a [X? {FP a A v a ) - x v {F p a a AP)\ 

+F»»Al-F?>A^ (121) 

where 



We can show that: 



d p M^ u P = dpM^P = 



+Fr A» a ~F»PAP a (122) 
(123) 



The spin of the quark can now be expressed by the fol- 
lowing matrix element: 



1 



(S k ) aP = -etjk / d 3 r (B a \ M ij0 \B a ) 



(124) 



The calculated expectation value should equal 1/2 in 
units H. 

We now express the conserved quantities in terms of 
the gluon profile functions and the quark wave functions 
or source functions. From Eq. (|119p we get for the total 
internal energy of the quark system: 

ro 

E mt = J d 3 r (B a \ 9 00 \B a ) =E-^E+^- J drl(r) 

o 

(125) 

where ro is the radius of the quark system and J(r) is 
given by: 



J(r) = S F 3 - 5iF + Ji 



[(*)' 




1 

9 




~ 7* 





F A + ^FoK 



FiFs — F0F4 



+4 ( rF 4 + 3Fl K ) 2 " ( rF i + ZFiKf (126) 
After reduction this expression simplifies: 



J(r) = S a F 3 - 5iF + -r 2 





— ( F1F3 - f fa 



F? -F[ 2 . 



(127) 



The internal energy of the system is now expressed com- 
pletely in terms of reduced profile functions and therefore 



is finite even if the original profile functions are singu- 
lar. For the 0(a°) bag solution with (3 = —2/3 we get 
E mt = E - (32/9)E < 0. In Section E] we will ad- 
dress the problem associated with the negative character 
of this internal energy. 

The total momentum of the system is easily calculated 
and equals zero, as expected: 



P' 



d 3 r(B a \0*°\B, 



0. 



(128) 



Since we are considering a system in isolation, it would 
not make sense to assign a momentum to it. Its momen- 
tum only becomes relevant when we embed it in ordinary 
QFT and the quark interacts with other particles. 
Let us now consider the spin. We obtain: 



(Si, ),, : = ((T,, ) „ ^ I - J drr (S F - S1F1) 



■ J y (rl% - F 3 + 2F if) (rFl + 2F l K) 


ro 

7^-/7- ( rF o ~ F + 2F * K ) ( rF l + 2F * K ) } t] 2 " 1 



32 

9a, 



The first term already provides the exact answer, hence 
the remaining terms should add up to zero. After reduc- 
tion we can convert this to the condition: 



1G 
27 



dr 



{{'*• 



32 
'9aT 



F 3 )F1 



or after integration by parts: 




drr S0-F4 - S1F1 



3 



(130) 



(131) 



To 0(a°) this equation is automatically satisfied by the 
solutions X(r) — Y(r) = (3Er, e.g. by the bag solution 
with /3 = —2/3. Its satisfaction for this bag solution can 
also be related to the identity So = So, which guarantees 
the validity of Eq. (| 130[) to lowest order. The condition 
(|13ip can assist in fixing the four free parameters defined 
in the perturbative scheme in Section fVl 

There is another conserved quantity for single quarks, 
namely color. The reason that this is not usually seen 
as an observable, is because quarks are usually not con- 
sidered in isolation (physically they only appear in color 
singlets). However, since we are considering single quark 
solutions to the field equations, we clearly have to exam- 
ine this property as well. Color is associated with the 
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invariance of QCD under local gauge transformations. 
The relevant infinitesimal transformations are (18(: 



and 



ip — > ip' — ip — i— a; • Xip 



A p — > = A M — eA M x w d^w. 



(132) 



(133) 



the color and spin condition to formulate the following 
boundary conditions: 




o UE-E-W' 

1 1/3 + EU 
UE-E- l V' 



r -> r . (141) 



The current corresponding to this transformation is given 
by: 

C = ^ip^Xip + ^F^ x A„ - x F^". (134) 



Since 

it follows immediately that: 



d^C = 0. 



(135) 



(136) 



Hence, we can define the conserved color density (/i — 0): 



-(F to xA„-A„xF°") fc 



\B 6 



(137) 



The first term in the right-hand side of Eq. (|137j) equals 
the bare color, so that we have to demand that the re- 
maining terms yield zero. Since this leads to a constraint 
on an integral, this condition is not easy to implement. A 
simpler condition arises if we use the identity, Eq. (|135p . 
again. This leads to the condition: 



(C° k ) =-(B a \ I d 3 xd„F^\B 
9s J 



1 

9s 



= — (B a \ I dSF° n \B p ) = (a\^ 



2 



Writing this out in profile functions, we obtain: 



3 ~ 



FqK 



«8 



(138) 



(139) 



where ro is the radius of the quark bag for the bag solu- 
tion. After reduction this condition can be written as: 



X 



E 



Y 



E- 1 



(140) 



Notice the similarity with the spin condition Eq. (|131l) . 
As for the spin case, this condition is automatically sat- 
isfied for the 0(a°) solutions Eq.(|S7)). We can combine 



Since either E or E is singular at r = ro , we need to treat 
these boundary conditions as limiting processes, rather 
than as point identities. For the bag solution we have: 



E : 



(x - x ) 2 



tan 2 xq 



(142) 



near the surface x = xq. After expanding A" in a Taylor 
expansion for r — > ro, we then arrive at three conditions, 
which allow us to fix the perturbative parameters v and 
uj. Numerically, we find: 



-0.052a s 



1.092a. 



(143) 



The third condition is more tricky as it involves an ad- 
justment of the radius of the quark bag. A rough cal- 
culation indicates that we have to adjust x from 2.043 
to 2.129, however, a more complete numerical study is 
required to get precise answers. 

In this section we considered the observables which 
play the role of boundary conditions in our formulation. 
Other observables, such as the magnetic moment, arise 
from the interactions with external fields, and must be 
calculated with standard perturbative methods. Clearly, 
our formulation has little to add to such calculations 
except perhaps providing physical cut-off parameters in 
view of the finite size of the dressed quarks. 



IX. THE ROLE OF GENERAL RELATIVITY 
AND THE DETERMINATION OF THE MASS OF 
THE LIGHT QUARKS 

As seen in the previous section, the total QFT energy 
of the dressed quark system is negative. Within the con- 
text of QFT there is only one obvious contribution which 
could compensate for this negative energy and that is the 
Higgs field. We expect that the two finite real Higgs so- 
lutions, whose classical approximations could be charac- 
terized by the values (j) — iVA'V'V would have such an 
effect and lead to an overall positive energy, which can be 
associated with the mass of the quarks. However, to en- 
able this treatment, one needs to know the exact nature 
of the Higgs Lagrangian at the fundamental level and 
preferably its relationship to the massive vector bosons. 
These Higgs terms would add an additional set of equa- 
tions and profile functions and therefore add considerably 
complexity. It therefore seems more natural to look ini- 
tially at the trivial Higgs sector with a zero Higgs field. 



18 



Strictly speaking, this "trivial" Higgs sector is not com- 
pletely zero, as the quark source term would induce a 
final Higgs field even if we choose the small Higgs solu- 
tion in the cubic field equation. However, we believe that 
it would be a fairly good approximation to approximate 
this lower branch by a zero Higgs field. The consequence 
of this assumption is that we need to find other solutions 
for the following two problems (1) The absence of a scale 
defining parameter; (2) The negative energy of the sys- 
tem. A system with negative energy and unconstrained 
scale would contract without bound in order to minimize 
its energy, unless its scale is fixed, which it is not in QCD. 
To prevent this collapse we consider the consequences of 
general relativity (GR). When the energy approaches the 
Planck scale, the effects of GR become important, and 
the large magnitude of the negative internal energy will 
halt the implosion. To treat these GR effects in a QFT 
context is non-trivial, and brings us in uncharted terri- 
tory. However, by treating these effects perturbatively, 
we avoid some of the tricky problems. Since this inclusion 
of the effects of GR, in combination with the vacuum en- 
ergy, leads to a spectacular agreement with experiment, 
we feel that there is good evidence for the validity of this 
perturbative approach. 

The effect of GR can be represented through the met- 
ric tensor, whose spatial component y^g modifies the 
spatial energy integral: 



d 3 x e ° -> / e °, 



(144) 



where the energy density is represented by the compo- 
nent O ° (depending on the metric convention, we may 
also need an additional minus sign). The spatial met- 
ric 3 g is controlled by the local energy density inside the 
quark. The simplest approximation is to replace this en- 
ergy distribution by an effective mass M e f f at the origin, 
with M e ff equalling the total internal energy —(23/9)E. 
By expanding the resulting integral in G we obtain a 
well-defined integral: 

f / 2M eff G\ 3/2 n 

£in*« d 3 xii + — f-\ e ° 

-) e °, (145) 

where 7 = 23/9. At the end of this section we will also 
discuss results where the point mass at the center is re- 
placed by an integrated density distribution based on the 
internal quark wave functions / and g. Carrying out the 
integral we obtain approximately: 

E mt ^-jE(l- 2*2^) , (146) 



8xp, 



where we replaced 1/r by its expectation value (1/r) = 
E/Sxo where 5 — 0.6019 • • • for the bag wave function. 
This expression reaches its minimum for 



E = 



' Sxq 
97G 



9x ~/G 



(147) 



The correction 3jGE 2 /5xq to the Minkowski metric 
equals 1/3, which means that the first order (pertur- 
bative) usage of GR is consistent. Both the radius ro 
and the internal frequency E are of Planck scale, with 
r = 8.8 lpi anck = 1.4 x 10 Ll9 fm and E = 0.23M Planck = 
2.8 x 10 21 MeV. The negativity of the internal quark 
energy was essential for the stabilization of the system. 
Hence, a property that appeared like a serious problem of 
the model actually was necessary to stabilize the system 
and to fix the scale. However, the relationship between 
the negative internal energy of Planck mass magnitude 
and the quark mass, which is positive and lies in the MeV 
range, is still unexplained. 

In order to address this problem we introduce the vac- 
uum energy density e, or what is equivalent a finite cos- 
mological constant. Recently, a cosmological model was 
developed which is based on the presence of this vac- 
uum (dark) energy |4|. The associated non-perturbative 
metric of the vacuum universe has a very simple time 
dependence: 



9nu = rinu9(t) = Viiv(ts/tf 



(148) 



where the characteristic time t s follows from the pos- 
tulated vacuum energy density e: t s — (S/SnGe) 1 / 2 . 
Another consequence of this cosmological model is that 
t 3 = I/Hq, where Hq is the Hubble constant, while t s also 
represents the age of the universe as measured by a co- 
moving observer. Fitting recent supernovae data one ar- 
rives at a vacuum energy density of e = 3.97 x 10~ 47 GeV 4 
and an age of the universe of t s = 13.8 billion years. The 
vacuum metric leads to the effective time-dependent vac- 
uum energy density: 



(t) = VWe = (t s /t) 3 



(149) 



We now make the novel assumption that the creation of 
a quark is associated with the creation of a small vacuum 
universe (a small big bang). The creation of this bag ter- 
minates at a time t c , when the size of the bag matches the 
size ro of the QFT quark bag. Since this creation time t c 
is very small, the effective vacuum energy density is very 
large (just like it was t c seconds after the big bang) , and 
the resulting energy can potentially cancel the internal 
energy, leaving a net quark mass of much smaller magni- 
tude. According to the Heisenberg uncertainty relation 
[24| the total energy during this time interval t c is un- 
certain by an amount h/2t c , so that it is possible during 
this time interval to create a mass of magnitude h/2t c . 
Hence, we set: 



h 1 
2.Z 



(150) 



and verify whether this assumption for the quark mass is 
consistent with the energy balance: 



m„ 



h 1 
2.Z 



(151) 
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where Ei nt is the minimum energy determined earlier. 
Since the quark mass will turn out to be infinitesimal 
compared to the other terms, it is natural to ignore the 
left-hand term: 

47T / Xq \ ^ fts^\ 
Eyacuum + E{ nt — — ^ — J € I — J + Ei nt ~ 0. (152) 

Clearly, the required cancelation between internal and 
vacuum energy can only happen if the internal QFT en- 
ergy is negative, so the negativity of the internal energy 
is again essential for the construction of a realistic quark 
model. The cancelation condition leads to a creation time 

to- 

This expresses the creation time - and therefore the mass 
of the quark - exclusively in terms of the cosmological pa- 
rameters G and e (or Hq), as the remaining parameters 
(x = 2.0478 • • • , 7 = 23/9 and S « .60) are fixed by the 
formalism. We obtain t c = 1.04 x 10~ 22 s. According to 
Eq. (1150)) this corresponds to a mass of the dressed quark 
of m q = 3.17 MeV, confirming the smallness of this mass 
in comparison to the Planck scale. Expanding the inte- 
gral in Eq. (1145[) to 0(G 2 ) only gives a small correction: 
one finds m q = 3.50 MeV (since there are other terms 
contributing in 0(G 2 ) which are not examined, we do 
not claim that this result is more accurate than the first 
order result). Replacing the mass term Mjr in the metric 
factor by a continuous distribution based on the internal 
quark wave functions one finds m q = 4.05 MeV. These 
three results give a good indication of the uncertainty in 
the result m q — 3.17 MeV, due to the uncertainties in 
unifying QCD and GR. Let us now compare this result 
to "experiment". 

Reference [25[ quotes the following lattice values for 
the average up and down mass: 3.8 ± 0.8 MeV, while the 
average value excluding the lattice is quoted as: 4.4 ± 1.5 
MeV. Very recently Davies et. al.[2l[ reported on a new 
determination of the light quark masses of much higher 
precision. They quote the following average values for the 
light quark masses: m q (2 GeV) = 3.40 ± .07 MeV and 
m q (3 GeV) = 3.07±.06 MeV, where the energy in brack- 
ets is the scale used in their calculation. The agreement 
between our result 3.17 MeV and these "experimental" 
values is extremely good, especially if one takes into con- 
sideration the enormous differences between the particle 
physics scale and the cosmological scale, from which our 
masses have been derived. Naturally, additional effects, 
like the electro- weak forces (which splits the degeneracy 
between u and d-quarks) and 0(a s )-terms (if a s ^ 0), 
have to be included in more detailed calculations. How- 
ever, the current level of agreement with experiment is 
clearly a strong endorsement of the theory presented here 
and confirms that this solution corresponds to the first 
generation of quarks. 



X. SUMMARY AND DISCUSSION 

In this paper we formulate a new set of methods in 
QFT that can be used to describe the internal proper- 
ties and masses of elementary particles. The exact non- 
perturbative operator structure of the fermion and bo- 
son quantum fields is constructed by demanding the sat- 
isfaction of the full set of coupled quantum field equa- 
tions. Starting from the basic QCD Lagrangian with 
bare pointlike massless quarks, we derive the internal 
wave functions of the dressed quarks and the form of 
the binding potentials. The binding mechanism of the 
localized quark-gluon state is topological in nature as it 
enabled by the non-linear structure of the equations and 
is independent of the magnitude of the strong coupling 
constant. Since all equations and fields refer to a single 
space-time coordinate, an additional operator prescrip- 
tion is required to reflect the opposite ordering of parti- 
cles and anti-particles. This prescription (the so-called 
R-product) also resolves the so-called cosmological con- 
stant problem, as it leads to a more fundamental and less 
ambiguous definition of the vacuum in QFT. 

The calculation of the quark mass is complicated by 
the fact that the QFT expectation value of the energy 
operator is negative. It is possible that by including an 
elementary Higgs-quark interaction, this problem may be 
resolved, however, this can only account for two of the 
three quark generations. The third - and probably basic 
- generation would then correspond to the trivial (zero) 
Higgs solution. However, in the absence of a Higgs field 
we need another entity to counter the negative QFT en- 
ergy and give the system a scale. We propose a novel 
solution for this puzzle by involving the perturbative use 
of gravitational forces at the Planck scale. In this way 
we arrive at a positive mass which is amazingly close to 
accepted phenomenological values for the light quarks. 

The current theory would get further support if we 
can deduce the existence of the two higher generations of 
quarks. However, as we indicated above, this may well 
require the inclusion of a Higgs quark interaction, which 
will entail more profile functions and a more involved 
formulation. But it also offers a chance to verify various 
models of this force and the associated electro-weak in- 
teraction at a more fundamental level than the standard 
model, by comparing the predictions to experiment. Af- 
firmation of the proposed Lagrangians would undoubt- 
edly lead to a more unified picture of Nature. Another, 
desirable extension is the perturbative inclusion of QED 
to explain the mass differences between quarks within one 
generation. Clearly, QCD cannot describe the dressing of 
leptons, so that for this sector the non-linear electro- weak 
interaction must induce the binding mechanism. The ap- 
plication of our formulation to unbroken SU(2) is fairly 
straightforward (for example we already found that the 
factor 32/9 is replaced by 3/2), but the chiral nature of 
the electro- weak interaction is more difficult to include. 
Again, this application to leptons offers a good opportu- 
nity to test various models of the electro-weak interac- 
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tions at the fundamental level, as it is possible that the 
complicated multiplet structure in the standard model is 
absent at a more fundamental level. The smallness of the 
neutrino masses would put extremely strong constraints 
on the possible theories, and if a model was successful 
in predicting the lepton masses then it would mark a 
historic breakthrough in our understanding of Nature. 

The application of the new methodology is not com- 
pletely limited to the description of the internal dynam- 
ics of elementary particles. It can also be used to con- 
struct free boson propagators using the new quantization 
method for solving the (homogeneous) equations of mo- 
tion. The R-product also could play a role in scattering 
diagrams if there are non-linear vertices with common 
space-time coordinates (as pointed out repeatedly, the 
R-product does not affect the algebra of operators de- 
fined for different space-time coordinates). In scattering 
calculations the Unite size of dressed fermions (emerging 
from the current theory) could also play a role in defin- 
ing physical cut-offs, and thereby ensure the finiteness 
of otherwise infinite diagrams. However, the size of the 
light quarks (about 8 Planck lengths) is beyond measure- 
ment. If the higher generations of quarks are indeed due 
to coupling to an elementary Higgs held, then their sizes 
would be substantially larger, and of the order of 0.004 
fm on average, so they might be measurable. Outside 
of these specific effects, the standard scattering formu- 
lation would remain unaltered, and so all the successful 
scattering applications of the standard model seem fully 
consistent with the current generalizations of QFT. 

Appendix A: The R-product 

The formulation of QFT contains an inherent asymme- 
try between particles and ant i- particles, which reflects 
the inadequacy of the mathematical language in which 
we express this formulation. Since we are constrained by 
this language we have to add another rule to QFT, which 
together with the usual rules constitutes a symmetric rep- 
resentation of the world of particles and anti-particles. In 
this appendix we define this rule. It should form an in- 
herent part of the foundations of QFT. Its absence in 
the current practice of QFT has led to many misconcep- 
tions and is also responsible for the so-called cosmological 
constant problem, which constitutes the biggest discrep- 
ancy ever between theory and experiment in the history 
of physics. 

Consider a typical operator expression in QFT-.TpOip. 
If we evaluate the vacuum matrix element of this operator 
and concentrate on the particle contributions we get: 

< 0|$Cty|0 >= ]T < 0\b%\0 > (4> a O<P p ) = 0, (Al) 

as expected. If we look at the expectation value for a 
particle state we get: 

< 0\b^Oi>\b\\0 >= (0 7 O^ 7 ), (A2) 



expressing the desired link between the state considered 
and the matrix element. However, when we consider the 
anti-particle contribution of the operator we have: 

< o|w|o >= < 0K4l° > (foot}) 

a,P 

= £(&0#), (A3) 

a 

which is in general non-zero, and can even be infinite, as 
the sum is unrestricted. Clearly, this asymmetry is unac- 
ceptable and illustrates the incompleteness of the stan- 
dard QFT description (notice that there is no problem 
with mixed terms such as b^dt, which are responsible for 
the quantum fluctuations). The origin of this asymmetry 
between particles and anti-particles is that the right-hand 
operator ip contains particle annihilation operators and 
anti-particle creation operators. This particular combi- 
nation is required by baryon (or lepton) number conser- 
vation. However, it leads to the possibility of creating 
anti-particles - and the impossibility to create particles 
- when operating with ip on the ket-vector, even if it is 
vacuum state. Similarly, the operator on the left, -0, con- 
tains a combination of particle creation and anti-particle 
annihilation operators, again required by baryon number 
conservation. This combination displays a similar un- 
acceptable asymmetry when operating to the left on a 
bra- vector. Let us see how we can restore the symmetry 
between particles and anti-particles. 

To eliminate vacuum terms, like those in Eq. (|A3[) . one 
usually applies the normal product to the bilinear expres- 
sion: 

d a dp -K dad^ := -dj 3 d a . (A4) 

Clearly, the anti-particle contributions to the vacuum 
matrix element now also vanish. The usual justification 
for this heuristic procedure in QFT is that in QFT one is 
not interested in absolute energies, but rather in relative 
energies, such as excitation energies. By consistently re- 
placing bilinear expressions by the normal product, one 
then eliminates such sources of infinities from the start. 
However, within GR energy leads to a modification of 
the metric, so that one cannot arbitrarily rescale the en- 
ergy. Thus the heuristic replacement of bilinear expres- 
sions by their normal product is no longer seen as ac- 
ceptable. This has led to the consideration of QFT vac- 
uum contributions to the energy as a real phenomenon 
(especially in cosmology), although this is in total dis- 
agreement with experimental observations. As we saw 
above, the normal product is needed to restore the sym- 
metry between particles and anti-particles, so we cannot 
simply dismiss it because it is a heuristic trick. Rather, 
we have to explain the origin of this heuristic prescrip- 
tion and should definitely not accept the reality of the 
vacuum contributions as a fait accompli, especially since 
these contributions are associated with an unacceptable 
asymmetry between particles and anti-particles and are 
in strong disagreement with observation. 
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The key towards the solution of this problem is to ex- 
amine the normal product from a different perspective. 
The normal product is defined as a rearrangement of the 
operators, so that the creation operators are shifted to 
the left of the annihilation operators. The normal prod- 
uct is also widely used in the derivations of the scattering 
expansion in QFT, although there it is used as a conve- 
nient mathematical construct and not employed heuris- 
tically. In its current heuristic application, Eq. (|A4l) . it 
is equivalent to the inversion of the order of the anti- 
particle operators, as the normal product leaves the par- 
ticle operators unaffected in the QFT expression ipOip as 
: bJJ)p := b^bp. Hence, an equivalent definition of the 
normal product for the relevant QFT expressions is that 
the order of the anti-particle operators must be reversed. 
We will refer to this re-ordering as the R- product. We 
will argue that there is a fundamental theoretical justifi- 
cation for this product, which is lacking in the heuris- 
tic application of the normal product. Hence, in the 
bilinear expressions considered one can justify the use 
of the normal product indirectly via the justification for 
the R-product. For longer expressions the usage of the 
R-product, rather than the normal product, leads to dif- 
ferent results. However, most of the long expressions in- 
volve different space-time coordinates, where neither the 
R-product, nor the heuristic normal product applies. If 
we analyze the usage of the heuristic application of the 
normal product in the existing QFT literature, we see 
that it is customary to define the (interaction) Hamilto- 
nian heuristically as normal ordered (see e.g. [26|), be- 
fore it is inserted in the matrix elements of the S-matrix 
expansion. This means that there is no heuristic normal 
ordering between operators corresponding to different co- 
ordinates x and x' (it makes no difference if we let x — > x' , 
for the current practice it is only important that the co- 
ordinates were different originally). In fact the heuristic 
usage of normal ordering between different spatial inte- 
grals in the S'-matrix expansion would destroy the whole 
structure and meaning of the S'-matrix expansion. Hence, 
the heuristic use of the normal product must be limited 
to strings of operators that all belong to a single coor- 
dinate x (we will later discuss the fundamental origin of 
this rule). But for such short expressions, correspond- 
ing to a single x, the normal product and the M-product 
are often equivalent, so that either one can be used. This 
equivalence for the simplest QFT expressions is the likely 
reason why the R-product has not been discovered before 
in QFT, as the usage of the normal product was often ad- 
equate. 

For longer expressions belonging to a single space-time 
coordinate the usage of the R-product leads to very differ- 
ent results than the normal product, so here it becomes of 
importance to use the M-product. In the current bound- 
state application of QFT we encounter many long chains 
of creation and annihilation operators belonging to a sin- 
gle coordinate x, so our transition to the R-product has 
major consequences for the current formulation (in stan- 
dard scattering applications of non- Abelian QFT it might 



also have consequences, as these also involve strings of op- 
erators at the same space-time point). Consider a longer 
local QFT expression, like dabpb^^ and b^dgdabp^t is 
essential to use expressions that actually occur in the us- 
age of QFT, since we are trying to correct a defect of 
QFT that occurs in its actual representation). Using the 
normal product we get: 



whereas the 



: d a bpb\d\ := b^dabp, 
: b\d\d a bp := b\d\dabp, 
-product leads to: 

d a bpb\dl = -bpb^dlda, 



(A5) 



b^dldabp 



-b\bpd a d\. 



(A6) 



Under the normal product the distinction between the 
two initial expressions is lost, while the R-product main- 
tains the integrity of the expression, while still ensuring 
that the vacuum matrix element vanishes. If we apply 
the resulting expressions to a one-particle state then we 
get zero under the normal product, while the M-product 
gives a finite outcome. Hence, the R-product leads to the 
survival of non-linear terms in the field equations which 
are essential for the existence of bound states, such as 
considered in this paper. 

The physical justification of the R-product is that it 
embodies Feynman's interpretation of ant i- particles as 
particles moving backward in time. In the S-matrix 
expansion the ordering of operators defined at different 
times is taken care off by the time-ordering procedures of 
standard QFT. However, when anti-particles all belong 
to the same time t, this ordering cannot be implemented 
by explicit time ordering. That is why one has to impose 
the R-product on local expressions to ensure the correct 
ordering of interactions and vertices. The terms forward 
and backward in time may seem a bit inappropriate for 
local operators defined at a single t, so let us analyze this 
concept in the current context. In the matrix elements 
the ket-vector is usually seen as the initial state, while 
the bra- vector is seen as the final state. From this per- 
spective, the operator right-most in a chain first operates 
on the ket-vector, and is "earlier" in time than operators 
further to the left. For a chain of anti-particle operator 
the left-most operator must first act on the initial state 
(i.e. the ket vector), etc. This can be accomplished by 
applying the R-product to expressions dependent on one 
space-time coordinate. 

One consequence of the R-product is that one cannot 
insert a complete set of states somewhere inside the oper- 
ator expression. Hence, the manipulation and reduction 
of an R-product requires special care, as we will also see 
in the following example. Consider the following reduc- 
tion: 



'^dJpd^bs 



= -b^d^d^bg = 



hp 



. (A7) 
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We see that if we commute the anti-particle operators un- 
der the R-product, we get an extra minus sign in the anti- 
commutator, compared to the original anti-commutator, 
Eq.(|lip. This is conveniently summarized by the follow- 
ing anti-commutation rule: 



[■■■{dldp} 



(-5 a p)---],x = x', (A8) 



where x and x' are the coordinates belonging to the two 
anti-particle operators. Since, all operators have to be 
expressed in R-product form, we have to apply this extra 
minus sign whenever we contract anti-particle terms cor- 
responding to the same space-time point in the equations 
of motion or in the expressions for the observables. These 
minus signs are absolutely necessary for deriving consis- 
tent equations, in fact we established the need for these 
signs, before we had discovered their origin in the R- 
product. We will see in Appendix IB) that the R-product 
also enables us to derive the general operator solution of 
the equations of motion by including operator chains of 
any length in the expressions. 



Appendix B: General operator solution of the 
equations of motion 

Since we are only interested in the creation and an- 
nihilation operator aspects in the current appendix, we 
will use a simplified, partly symbolic, notation. Consider 
first the Dirac equation. Writing this in the form: 



H i> = i/>A, 



(Bl) 



we get a series of one- and many-body terms on the right, 
when we use the expansions Eq. ((SJ and (|9]) for the quark 
and gluon field operators: 



E 



(B2) 



After contraction, some of the right-hand terms lead back 
to the single particle annihilation or anti-particle creation 
operator. For example, 



b a bpb 7 



5 a pb 1 \ 



dLdab. 



-Safjb^, 



(B3) 



where the minus sign arises from the R-product. How- 
ever, there are also a lot of many-body operators. For 
example, we have: 



bp b a b.y ; 



didpb 7 



-dpdlfi y , 



(B4) 



both of which correspond to many-body terms. Notice 
that when reducing operator expressions in R-product 
form, the anti-particle creation operators have to end up 
towards the right after contraction. After, carrying out 
the R-product in the final reduced expression, these cre- 
ation operators will then end up to the left, as desired in 



reduction processes. This explains why the second term 
in Eq. (|B4p must be considered a many-body term. 

Writing out Eq. (|B2j) explicitly, we can cast the quark 
field operator in the following form: 



1> 



J2 Ma + £ 4C + E [ b i b ^€ 

a, /3,7 



pap 



(B5) 



The reason that we have 6 many-body terms instead of 
the expected 4 x 2 = 8, is that two terms, namely b^bpdl 
and b a dpdL, receive contributions from two products. In 
order to have unique profile functions we demand that 



« = -o and c; 7 = -^x- 

(B6) 

For the gluon field operator we can carry out a similar 
expansion, by demanding that it has the same form as the 
many-body source function. For the quartic component 
of the gluon field operator we now get 16 terms, minus 
the duplicate products, i.e. 12 in total. Again, we have 
to impose anti-symmetrization conditions on the gluon 
profile functions, to ensure their uniqueness. 

To simplify matters we will concentrate on the particle 
terms, suppressing the particle indices for the moment. 
We then have: 

A = E A <*P b l b P + E A <»Pi5 b l b \ b i b 5, (B7) 



a/3 

where we demand that 



A 



-A 



-A 



a/3 #7 



.4 



(B8) 



The Dirac equation for the b a bpb 7 component then reads: 
Ho4>af3~ ( = —~<l>pA ai + -4>~fA a p + 2(f) t A ta p 1 

+(j>aPeA e ~/ — 4> al<L A el + 2(f) aeT A T(i p~ / . (B9) 

Hence, even if we limit ourselves to particle terms, we still 
get fairly complicated equations. Because of the first two 
terms on the right we do not have a solution </> Q ^ 7 = and 
A a /3-yS = 0. Hence, the solution of the field equations nec- 
essarily contains many-body components. Iterating the 
equations further, we can use the same argument in the 
next order, and so on, ad infinitum. Hence, the solutions 
of the field equations are operators that contain compo- 
nents that are of arbitrary high order in the creation and 
annihilation operators. What is more, these components 
cannot be treated perturbatively, as they are the same 
order as the one-body terms. This series is cut off in 
practice by the finite size of the wave function space, as 
products containing two identical creation or annihila- 
tion operators gives zero. For scattering problems the 
wave function space is infinite, however, here one often 
deals with operators at different space-time points be- 
tween which the R-product does not feature. 
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After analyzing Eq. (|B9[) and the corresponding quartic 
gluon operator equation, we find the following solution: 



(BIO) 



and 



(Bll) 

We can use Eq. (|B10|) to write the quark field operator as 
follows: 



a/37 



S K ^ = ( 1 - N )J2 b ^ ( B12 ) 



where 



N = Y,b%, 



(B13) 



is the particle number operator. Similarly, we can write: 
i = ^6l(l-%V (B14) 

a/9 

We can show that the non-linear terms in the gluon field 
equations do not invalidate this result, by noticing that: 

AA ->■ &t (1 _ N)bfiA a pb\ (1 - iV^A^ -► 

& a (l-iV)Ma^5 + --- (B15) 

where the • • • stand for higher-order terms, which can 
be ignored at this stage. So products of A's also lead to 
terms of the form, Eq. (|B14[) . Their functional form is the 
same as the one-body terms considered, which ensures 
that the scalar field equations have a simple connected 
structure without need for perturbative solutions at the 
operator level. 

We can extend this analysis to anti-particles. This 
analysis is a bit trickier because of the implied R-product. 
We find: 



^ = (i-Jv-w)^(w„ + 4< 



(B16) 



and 



A = 



£ [&t (1 — N — N)bpA™ + d a (l -N- N)b A a £ 

a/3 

+&t (1 -N - N)d\A^ + d a (l -N- N)d\A a ^ 



where 



N 



Y^ded\, 



(B17) 
(B18) 



is the effective anti-particle number operator under the 
R-product. The simplicity of this result is only true under 
the R-product. Once, we carry out this product explic- 
itly, the elegance and simplicity is lost. 

We can now consider terms of higher order, some 
of which already appeared in our expressions, e.g. in 
Eq. (IB15[) . Surprisingly, we can get an exact solution of 
the many-body equations to any order. We can write: 



a 

Here we defined: 



(B19) 



(1-N-N) 



2- N - N n-N - N 



72- N -N 



(B20) 



where we generalized the definition of the factorial to 
allow for operators. The gluon field operator solution 
becomes: 



a/3 



ap 



(B21) 



where we re-introduced the upper and lower indices in 
A. The amazing thing is that we have found an exact 
solution without increasing the number of profile func- 
tions. These general solutions of the operator problem 
allow us to reduce the operator equations to scalar equa- 
tions. The operator T,^ is a vacuum projection operator, 
as it yields zero unless N = N = 0. However, within op- 
erator expressions it is surrounded by single creation and 
annihilation operators, which implies that in practice it 
acts like a one-body projection operator. Hence, we have 
proved that this exact solution of the operator field equa- 
tions only yields solutions when operating on one-body 
state vectors. Furthermore, the quark creation operator 
operating on the vacuum, which originally created a bare 
quark, now creates a dressed quark state, and not just a 
component thereof. It would seem natural to represent 
many-body components (like sea quarks) in the state vec- 
tor. However, as this study has shown, these many-body 
components appear in the field operators, and conspire 
to leave the state vector simple, so that a (dressed) quark 
still can be represented by a single creation operator and 
the simple state vector b' a \0). 

Another consequence of this result is that these self- 
consistent solutions of the field equations cannot describe 
multi-quark states, such as the proton or the pion. This 
is a bit disappointing, since the emergence of the MIT- 
bag mechanism from QFT, would suggest that a similar 
derivation can be used to justify existing phenomenolog- 
ical approaches to the multi-quark proton case. How- 
ever, within the present operator solution this is only 
feasible if the multi-quark state is treated as a one-body 
state (as the proton was before quarks were discovered) . 
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But an "elementary" proton is a color singlet, and thus References 

it would no longer be possible to treat it via the non- 

Abelian SU (3) -theory, as was essential for the emergence 

of the bag solution. So the existing techniques, such as 

lattice gauge calculations, remain the tool of choice to 

deal with multi-quark systems. 
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